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INFLUENCE MATRIX PROGRAM FOR 


AERODYNAMIC LIFTING SURFACE THEORY 


Richard T, Medan 
Ames Research Center 
and 

K. Susan Ray 

Computer Sciences Corporation 


ABSTRACT 


This document is a description of and users manual for 
a USA FORTRAN IV computer program which computes an 
aerodynamic influence matrix and which Is^ one of 
several computer programs used to analyze lifting/ thin 
v/ings in steady, subsonic flow according to a kernel 
function method lifting surface theory. The most 
significant features of the program arc that it can 
treat unsymmet r ica 1 wings, control points can be placed 
on the leading and/or trailing edges, and a stable, 
efficient algorithm is used to compute the influence 
matrix. 
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INFLUENCE MATRIX PROGRAM FOR 


AERODYNAMIC LIFTING SURFACE THEORY 


Richard T, Medan 
Ames Research Center 
and 

K, Susan Ray 

Computer Sciences Corporation 


1 INTRODUCTION 


This document is a description of and users manual for 
a USA FORTRAN IV computer program which computes an 
aerodynamic influence matrix and which is one of 
several computer programs used to analyze lifting, thin 
wings in steady, subsonic flow according to a kernel 
function method lifting surface theory. The most 
significant features of the program are that it can 
treat unsymmet r i ca 1 wings, control points can be placed 
on the leading and/or trailing edges, and a stable, 
efficient algorithm is used to compute the influence 
ma t r i X , 


Questions concern inc either this document or the 
computer program or the associated computer programs 
should be directed to; 

R. T, Medan 
Ma i 1 Stop 247-1 
Ames Research Center 
Moffett Field, 

Ca. 94035 
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2 PROGRAM DESCRIPTION 

This program implements the algorithms discussed in 
ref. 1 for determining the aerodynamic influence matrix 
of a thin> lifting, planar wing in linearized, steady, 
subsonic flow. (i.e., the downv/ash modes induced by 
/he pressure modes). The program utilizes the new 
technique (discussed in ref. 1) which requires fewer 
quadrature points yet still calculates the influence 
functions accurately enough to guarantee convergence 
with an increasing number of spanwise quadrature 
points. 

The types of wings which may be treated by the program 
include asymmetric ones and ones with curved edges. 
Also control points on the wing leading and/or trailing 
edges may be handled. The v;ing geometry and spanwise 
integration points are input to the program from a disk 
file or tape (see section 5.1) created by the geometry 
program (ref. 2). The influence matrix vyhich this 
program calculates is written on a disk file or tape 
(AIM file, see section 5.2 for a complete description). 
This AIM file is then used by the equation solving 
program (ref, 3), which determines the coefficients in 
the expansion for the lifting pressure coefficient, 
Delta-Cp. Also the AIM file may be used by the 
boundary condition program (ref. U) and the planform 
plotting program (ref. 5) since the wing control points 
may be changed by this influence matrix program. 

This program was designed to operate under a command 
format. After some initial input is entered, the 
program prints a + sign (in the conversational mode). 
At this point the user enters a command. Then the 
program performs an action associated with the command. 
This action may consist of reading a variable or 
variables, reading In the geometry of a new wing, or 
beginning the computation of the AIM file for the 
current wing. 

This command format approach has been found to be 
extremely flexible and to make the program very easy to 
use. For example, it allows the user to enter data in 
any order that he wishes and allows a conversational 
user to correct mistakes easily since each command may 
be given as many times as the user wishes, A complete 
description of the available commands is given in 
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section 4. 

Another particularly useful feature of this program is 
that, with a single program run, the user can create 
nxjre than one AIM file. These files can be AIM files 
for different v/ings and/or AIM files for the sane wing. 
AIM files for the same wing can differ from one another 
because of changes in the Mach number or else different 
choices for the number of spanwise integration points 
or other integration parameters or else because of a 
different set of control points. The ability of the 
program to change the chordwise control points 
arbitrarily and the spanwise control points to a 
limited extent from those (default) control points 
stored on the geometry file has been found to be 
another very useful feature. 
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3 USERS INSTRUCTIONS 

3.1 INITIAL SETUP AMES' TSS SYSTEM 

For either batch or conversational processing the 
following TSS commands must be given. These commands 
are required once and only once for each user ID, The 
first three commands create the identification number 
file named IDFILE. This file contains four zeroes in 
binary form, 

SHARE MEDAN, FSARTM, INIDFILE 
COS MEDAN, IDFILE 
DELETE MEDAN 

SHARE MEDAN, FSARTM,LSPROG, VI 

3.2 CONVERSATIONAL USE ON AMES' TSS SYSTEM 

All integer data should be entered in a 1615 format and 
all floating point data in 8F10.0 format, 

USER; After logging on enter the following; 

AMES USYSLIB 
JOBLIBS SYSULin 
JBLB MEDAN 


It is not necessary to issue DDEFS for the 
geometry, identification number, and influence 
matrix files because the program issues its own 
DDEF commands to TSS, 


USER; 

CALL INF$ 



PROG; 

ENTER BATCH 



USER; 

Enter carriage return 

for conversational mode. 

PROG: 

ENTER IDl. 



USER: 

Enter identification 
file. If a negative 
most recent geometry 

number 
numbe r i s 
file w I 

of the geometry 
entered, then the 
11 be used. Enter 


zero to terminate execution. 


PROG; + 
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USER: At this point the user must begin entering 

commands. After each command he should enter 
the Input, if any, associated with that command. 
The program does not prompt for this input. If 
the command Is any of the commands except STOP, 
START, or WING, the program will respond with 
another + sign and the user must enter the next 
command. With a STOP command, the program will 
terminate, IVith a START command the program 
will request a value for ODISK, If output on 
the terminal is desired, enter zero. Otherwise 
enter an integer from 1-9. The program will 
respond by giving the name of the file where the 
output will be found. Then the program will 
commence with the calculation of the influence 
matrix. This may take many minutes on the AHps 
IBM 360/67 since the calculation can be lengthy. 
After the calculation has been finished, the 
program will respond with another + sign. At 
this time the user may make mod If 1 cat ions to the 
data and calculate another influence matrix or 
else he may give the STOP command. After a WING 
command the program loops back to the point 
where it requests I Dl , After a WING command has 
been given, the program resets variables which 
have default values. This moans that, for 
example, if the user had given the NSTOPF 
command once, he would have to give it again 
after a WING command to retain the effect of the 
NSTORE command, 

3,3 TSS BATCH MODE 

The batch mode operates the same as the conversational 
mode v\/ith the sole exception that a "T" must be put In 
column 1 on the first card (for the logical variable 
BATCH). This "T" suppresses all subsequent 
conversational prompts. 
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4 DESCRIPTION 0£ COMMANDS 

In all cases the first three letters of a commnnri are 
sufficient input. All integer data should be entered 
in a 1615 format and all floating point data in 8F10.0 
format. The inputs if any^ associated with each 
command is to be entered on the following line in 
conversational processing and in batch processing is to 
be on cards immediately following the command card. 
The same command may occur more than once. This is 
useful in correcting data entered in error. 

Although the commands nay, for the most part/ be 
ente red in any orde r the foil owing command sequences 
must be observed: PP , if it is to be given, must 

precede the NMAX and CV/TYPE commands. The START 
command must precede the VMNG and STOP commands. The 
MM command must precede the SWTYPE, KK, and MREF 
commands , 

Certain commands reverse their effect each time the 
command Is used. In other words, if the subject 
commands are given twice, the effect will be the sane 
as if they were never given. Also the variables 
associated with these commands are reset after a WING 
command. This means that if, for example, the NOHEAD 
command were in effect, the user would have to give 
this command again following a WING command in order to 
retain the effect of the NOHEAD command. The commands 
which reverse their effect are: IWRITE, NSTORE, IWCCP, 

IWSEP, SUPPRESS, CCI, NOHEAD, and MINT. 

Besides resetting the variables associated with the 
reversible commands, the WING command causes all other 
variables to be reset to their default values. These 
default values either come from the geometry file or 
are set by this influence matrix program. 

The commands which are the most Important and basic 
are: JJ, which adjusts the number of 5 Danv^fise 

integration points; MM and SWTYPE, which adjust the 
spanwise control points; PP and SV^TYPE, which adjust 
the chordv^ise control points; START, which initiates 
the major part of the computation; IVING, which reads in 
geometry data for a new wing; and STOP, which 
terminates the execution. These and all of the other 
commands are explained fully below: 
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Q£1 

Effect; Causes chordwlse convergence information to be 
printed* For a typical case this command may 
cause an sddltional 5000 lines of output, so 
this command Is not generally used except when 
debugging. 


CONTINUE 

Effect: Causes processing to continue in batch mode 

when an erroneous command is encountered. 
Othervi/ise^ execution will terminate. In the 
conversational mode an invalid conrrtand message 
will be Issued and processing v/lll continue.’ 

CWTYPE 


Input: CWTYPE (integer variable) 

(CHICP(P),P=1,PP) (if CWTYPE<0). 

Effect: Causes program to read CWTYPE, an integer 

denoting the type of chordwlse control point 
distribution. (CHICP(P),P=1,PP) will then be 
computed^ if CV/TYPE >,0 or read if CWTYPE<0. 
(CMICP) is an array of PP elements giving the 
relative chordwise location of the control 
points (O.lCHICPa.). The control points 
which are computed are the Wagner points if 

Multhopp points if 
CUTYPE>0 (fig, 2). The accuracy with which 
the downwash modes are calculated decreases 
for control points approaching but not on the 
leading and trailing edges. If this command 
is^ not given the default chordwlse control 
points on the geometry file will be used. 
This command must be preceded by a PP command 
If It is necessary to give a PP command. 


■ P^ m Q 

Input : 


DELTAO 


Effect: 


Causes program to read DELTAO, 
parameter in the spanwise 
described in ref. 1. It has a 


DELTAO is a 
in teg ra t i on 
default value 
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of AR**2 where AR is the aspect ratio, 

EPS 

Input: tPS 

Effect; Causes program to read EPS, an integration 
parameter discussed in ref. 1. The default 
value for EPS is .5. 

lV^iCC?_ 

Effect: Causes printing of chordwise control points. 

This printing will not occur until after the 
START command has been given. 

I WE1T£ 

Effect; Causes printing of the Influence matrix. 

i. i is c r. 

Effect; Causes printing of spanwise control points. 

This printing will not occur until after the 
START command has been given. 

M 

Input: JJ 

Effect: Causes a new JJ to be read. JJ is the number 

of spanwise integration points. The default 
value is JJMAX, which is a number read from 
the geometry file. If necessary, JJ will be 
adjusted downwards until M00( JJMAX+l, J J+1 ) =0 . 
This adjustment is necessary because the 
integration points used must be Identical to 
or a symmetric subset of the integration 
points which are on the geometry file. 


Input: KK 

Effect: Causes new KK to be read. KK is the maximum 

order of the spanwise pressure modes. The 
actual number of modes computed is equal to KK 
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um 

Input : 
Effect : 


LMIN 
Input : 
Effect ; 


MACH 
Input: 
Effect : 

m 


unless M0DES<.0 (see below). Th Is command must 
follow the MM command if KKH^M^ otherwise this 
command is not needed. 


LMAX 

Causes program to read LMAX, which is the 
maximum number of chordv/ise integration points 
which should be used (excluding the point on 
the leading edge). This number will be 
adjusted downwards until LMAX+l-( LMI N+l)*2** I 
where I is an integer. it will also be 
adjusted downwards until it is as small or 
smaller than LMXMX in subroutine INITAL. 
Currently LMXMX = 511. This relationship 

occurs because of the nesting of the chordwrse 
quadrature points (ref. 1). 


LMIN 

Causes program to read LMIN, which is the 
minimum number of chordwise integration points 
which should be used. If this command Is not 
given then LMIN will be computed according to 
the type of chordwise control points which 
will be used. If CWTYPE<0, then the program 
assumes LMIN=1, If CViTYPE>.0, then LMIN will 
be computed so as to align the chordwise 
control points with the chordwise integration 
points. The latter improves the convergence 
rate of the chordwise integration when the 
distance between the spanwlse control and 
Integration points is small. 


MACH (floating point) 

Causes program to read a new Mach number. The 
default value is the value stored on the 
geometry file. 
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I nput : 
Effect : 


Input : 
Effect; 


MREF 
Input : 
Effect ; 


NINT 
Effect : 

NMAX 
I nput : 


MM 

Causes program to read MM, which Is the number 
of spanwlse control points. The default value 
for MM is the value of NF stored on the 
geometry file. This command also sets KK 
equal to MM and MREF equal to MM, so. If the 
KK and/or MREF commands are to be given, they 
should follow this command. 


MODES 

Causes program to read MODES, which denotes 
the type of spanwlse pressure modes to be 
used. The default value Is 1, which causes 
both the symmetric and antisymmetric downwash 
modes to be calculated. If MODES<0 only the 
antisymmetric downwash modes will be 
calculated. If MODESTO only the symmetric 
downwash modes will be calculated. 


MREF 

Causes program to read MREF, which Is a 
reference number for the calculation of the 
spanwlse Integration points. MREF is the same 
as NN In the geometry program (ref. 2) and is 
explained therein. Also see the explanation 
of the SWTYPE command. If MREF=MM this 
command Is not needed. If MREF+MM then this 
command should be given after the MM command. 


Causes number of chordwise Integration points 
for each control point and spanwlse 
integration point to be printed out. 


NMAX 
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Effect: 


NQ H£AP 

Effect : 


NSTQRS 

Effect: 

PP 

Input : 
Effect : 

START 
Effect : 


ST^E 
Effect : 


Causes program to read NMAX^ which is the 
number of chordwise pressure modes. The 
default value for NMAX is PP. If NMAX+PP^ 
this command must follow the PP command. If 
NMAX=PP, this command Is not necessary. 


Suppresses the case information which is 
ordinarily printed after ODISK is entered and 
just prior to computing the Influence matrix. 
This command is generally only used^ in the 
conversational mode when the output is at the 
terminal and then it serves to decrease 
terminal time. 


Causes the program not to create an influence 
matrix file. This command is generally only 
used when testing the program and frequently 
is used with the I WRITE command. 


PP 

Causes program to read PP, which is the number 
of chordwise control stations. This command 
also sets NMAX equal to PP, so if the NMAX 
command is to be given, it should be given 
after the PP command. 


Causes the program to print a page of heading 
material containing the integration parameters 
and other information and then starts 
computation of Aerodynamic Influence Matrix. 
As the computation is being done the influence 
matrix is stored PP rows at a time (if the 
NSTORE command has not been given). 


Terminates execution. 
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SUPPREaa 
Effect : 

SVJTYPE 
I nput : 

Effect : 


I nput : 
Effect: 

WING 
Effect : 


Causes suppression of warning messages 
occurring because the chordwise integration 
fails to converge. 


SWTYPE (integer) 

(NINDEX(M),M=1,MM) (If SVJTYPE+O) 

Causes program to read SV/TYPE/ which is an 
integer denoting the type of spanwise control 
points. If SWTYPE=0, then the usual Multhopp 
distribution will be used. If SWTYPE+0, then 
the array (NINDEX) will be read and the 
spanwise control points will be given by 
COS(N INDEX (M) *3. 14159. . ./(MREF+D) for 

M=1^2^...MM, This command should follow the 
MM command, if the MM command should be given. 


A TSS command of 80 characters or less. 

The command is passed to the AMES* TSS 
operating system. After the system processes 
the command, control returns to the program. 
This command Is a special one for the AMES* 
TSS version of the program. 


Causes program to request a new value for IDl 
so that another wing can be processed. All 
variables and commands are reinitialized after 
this command is given* 
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5 SAMPLE TERMINAL SESSION 

A sample terminal session on the AMES' 560/67 T5S 
computer system is given in this section with 
additional comments added in parenthesis. During this 
session 3 aerodynamic influence matrices for an AR = 2 
v/ing were computed. 

For the first influence matrix the default control 
points and integration stations were changed. For the 
second influence matrix everything was the same as for 
the first except that the Mach number was changed. 
Both matrices were printed out, but not stored. 

For the third matrix the default integration stations 
and control points were used. The defaults were easily 
restored by using the VHNG command. This matrix was 
computed and stored/ but not printed because of its 
large size. 

The following example was reproduced from an actual 
terminal session: 

LOGON user id/password/ te rminal id 

AMES USYSLIB 

JOBLIBS SYSULIB 

JBLB MEDAN 

CALL INF$ 

(The influence matrix program is now in control.) 

ENTER BATCH 
(carriage return) 

ENTER IDl 
4 

+ 

IWRITE 

+ 

JJ 

47 

+ 

PP 

3 

+ . 

CWTYPE 

0 

+ 

MM 

3 
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+ 

IWCCP 

+ 

I WSCP 
+ 

NSTORE 

+ 

TSS 

CPUTIME? 

0.765 SECONDS 

+ 

START 

ENTER ODISK 
1 

OUTPUT IS ON .. .OUTPUT. AIM, Nl. . . 

ID2= 0 

(Now the program prints output and computes the 

Influence matrix.) 

+ 

TSS 

CPUTIME? 

2,978 seconds 

(By subtracting the previous CPUTIME from the one 
above, the user may determine how many CPU seconds It 
took to calculate and store the influence matrix 
program, ) 

+ 

MACH 

0.8 

•f 

START 

ENTER ODISK 
2 

OUTPUT IS ON ,. .OUTPUT. AIM. N2. . . 

ID2= 0 
+ 

TSS 

CPUTIME? 

4.944 seconds 

+ 

WING 

ENTER IDl 
4 

+ 

IWCCP 

+ 

U^SCP 
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TSS 

CPUTIME? 

5.355 seconds 

+ 

START 

ENTER ODISK 
3 

OUTPUT IS ON ...OUTPUT. AIM. N3... 

102= 5 

(The influence matrix will be on the fil 
AIM.X004.X005 .) 

TSS 

CPUTIME? 

59.450 SECONDS 

+ 

STOP 

TERMINATED; STOP 

(The operating system is now in control.) 

PRINT OUTPUT.AIM.N1^PRTSP=EDIT,STAT10N=RMT05 
PRINT BSN=????, ??? LINES 

PRINT OUTPUT.AIM.N2^PRTSP=EDIT^STATION=RMT05 
PRINT BSN=????^ ??? LINES 

PRINT OUTPUT.AIM.N3,PRTSP=EDIT,STAT10N=RMT05 
PRINT BSN=????^ ??? LINES 

LOGOFF 


e named 
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6 INPUT FILES 

The following disk files are read by the program. The 
Ames* TSS version of the program issues its own DOFF 
commands for the files/ so none need be given. For 
other systems appropriate control cards will have to be 
supplied for units 7 and 9. 

6.1 GEOMETRY FILE 

This file is a variable record length file and is read 
from unit 7, 

The first record contains identification and title 
Information Including the number of control points and 
integration points. 

The next record contains the chordwise control points, 
the array of indices from which the spanwise control 
points are derived, the tangents of the wing edge sweep 
angles at the integration stations, etc. For a 
complete description of this file see Ref. 2. 

On the AMES* TSS system this file has the name GEOM.XI 
where I is the numerical value of IDl. 

6.2 IDENTIFICATION NUMBER FILE 

This file is read from unit 9 and rewritten on unit 9 
and contains identification numbers in binary form. 
The second number on this file (ID2) Is incremented by 
1 and then the file Is rewritten using the Incremented 
value of ID2. I D2 is the identification number which 
will be used for the aerodynamic influence matrix file. 
It Is printed out with the heading information for each 
case . 

On the AMES* TSS system this file has the name IDFILE. 
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7 


INFLUENCE MATRIX (AIM) FILE 

This file is a variable record length file and is 
written on unit 11* The Ames* TSS version of the 
program issues its own DOEF command for this file, so 
none need be given. For another system appropriate 
control cards would have to be supplied for unit 11. 

The first record contains identification and title 
information plus information about the size of the 
matrix and location of spanwise and chordwise control 
points. 

The second and subsequent records contain the influence 
matrix itself. 


On the AMES* TSS system this file has the name 
AIM.XI.XJ where 1 is the numerical value of IDl and J 
is the numerical value of I D2 , 

The specific information stored on the file is the 
fol 1 owing : 

FIRST LOGIC AL RECO^ 


IDl 

1 dent i f i cat i on 

number 

of 

the 

geometry file. 

ID2 

1 dent i f icat Ion 

number 

of 

AIM 

file. 


NSEQ 

(TITLE) 

NTITL 

PP 

CWTYPE 

MM 

MREF 


This integer no longer serves a practical 
purpose . 

Titling information of up to 80 characters. 
This array is 26 integer words long. 

Number of words in TITLE to be printed. 

Number of chordwise control stations 
( integer) . 

An integer denoting the type of chordwise 
control point distribution. 

Number of spanwise control points. 

A reference number for the calculation of the 
spanwise integration points. 


17 



Influence Matrix Program 


SWTYPE 

NMAX 

KK 

MODES 

UNSYM 

LMIN 

LMAX 

JJMAX 

JJ 

MACH 

EPS 

DELTAO 

(CHICP) 

(NINDEX) 

(ETACP) 


An integer denoting the type of spanwise 
control points. 

Number of chordwise pressure modes. 

Maximum order of the spanwise pressure modes. 

Denotes the type of spanwise pressure modes 
whose downwash modes have been calculated and 
are stored on subsequent records of this file. 

Integer which/ when not equal to 0, indicates 
an unsymmetric wing. 

Minimum number of chordwise Integration points 
which was used. 

Maximum number of chordwise integration points 
which was used. 

Maximum number of spanwise integration 
stat I ons . 

Number of integration points which have been 
used , 

Mach number. 

An integration parameter. 

A parameter in the spanwise integration. 

An array containing PP chordwise stations at 
which the downwash modes have been determined. 

An array containing MM integers for computing 
the spanwise control point positions. 

MM spanwise stations at which the downwash 
modes have been determined. Although (ETACP) 
can be derived from (NINDEX) and MREF/ it was 
considered convenient to store these points 
despite the redundancy. 
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IIEXT MMP RECORDS 


( ((ALFA(rJ,K,P)^ri = l,NnAX),K = K.L,KK,KJUMP),P=l,PP) 

VJhere KL=1 if *iODES>.0 or KL = 2 if MnnF3<0 and KJI 1 MP = 1 if 
'■10DES>0 or KiinMP = 2 if MODESO and MMP='iM i^ the vnng is 
unsynr.iet r i ca 1 , or MMP = ( i-M+1 )/ 2 if the vn ng is 
symmetrical and ilODES>.0 or MMP = r'M/2 i^ t^e ‘.v/in<^ is 
symmetrical and flODES<0, 


ID 
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TITLF , s KECfftNr,ULA« Att s 2 

iDt = ^ . 

ID? s ft , 

mSEQ z . \ 

PP s . , 

CwTYPE s ft ■ 

MH s . 

MREF = . ? 

SWTYPE' aft 

NMAX = T 

KK a T ' 

mooes 3,1. 

l)N5Y^' =. \ ft 

LMIN a- .1 

I mAx s Ell 

JjM&y 3 i>J 1 

JJ =' ■ ^7 , 

bEiT AO = a , 00000 

F P 5 3 0 ,^0000 

MACH = .0,000 00. 

livF ST0‘ s , 1 . 

rcHir'p) 

O.OOCQOO 
n^^OftOftft 
i[ no 00 0 0 

J Ik:DEX t'fA 

1 I 0,70711 

2 P «0,00fti00 

3 3 -0,70711 


matrix . . 
11-13-73 



starting 2000 ^ Lonp 


rs3 


FT&= 0,70710^6 OTsO.OOOOOOO 

0‘.7?4i??63 -0 ,t^oQt7U 

CAI FA (N, p/ n , Nrl . NmAX) 

r,52lB7feO .0.^^66590 

( A L F A ( N , 3 , n f = U ^ ■ A V ^ 

l,33ai^S6 •n.lRa^SBX 

FTAs 0.70 71 0B6 T H J =0 , ii99909S 
f A I F A ( N' . 1 , 2 ) ^ N = 1 f A y 1 

r.l2lB96b 0, 2^00033 

fALFA(N, 2, ?) , I'JS 1 f ) 

2',5702220 0,707 1 1 1 2 

fALFAfN, 3, ?3,N«l,^^Ay^ 

2,556809^ 0,7^00026) 

FTA? 0.707 1 0B6 C HI si , 0000000 
C Al. F A ( K-, 1 , 35 , Nsl » 1 

r, 2 795 b 0 0 0,5909996 

CaIFACH, 2, 3),N=l,^^'Ay5 

3’,0‘5B933« 1 , 7 7 06^39 

fAt>A(*J. 3. 35 f X) 

3,10^191 1 1,6996695 

FTAe.O.COOOOOO Cr'T-0,OOOaOtiO 
c A I P A C N , 1 , n » S 1 » A X 5 

l',0?76797 -0.96293b3 

f A I F A f N . ? , n » ■" s 1 » ^ ‘ ' * X 5 

.0,0000000 0,0000000 
( A I F A c » 3 , n / = 1 * ^ ^ X ^ 

*1 ,^637335 0 ,3379765 


«Tnt,Pts,= 260 

"0t9i 97662 

-0,S35790l 
-0,3357690 
fiTNt,PTS,3 572 
o', 5995938 

0',9359369 
cr.5005l2l 
aTNT,PT5.= 260 
-0,91 9766 o 
-0,5357909 
-0.3357699 
dTMT,PTS,s 256 
-0,61 0 36 35 
0 0 0 0 0 0 0 0 
a',529lB22 



n^^s^.OpOC'ooooo 
fAI.F*’AC*^» 1* 2) » I p ) 

1,2415^^775 0,2500031 

(ALPA(N, a, 2) rNsUNMAy^ 

o'. 0000000 0,0000000 

(ALFA(N, 3, 2>,Nsl,fvMAy5 

.2,fe82(>989 .0,7500027 

( 

FT As .0.0000000 CHlsl.OOOOQOO 
(ALFA(f^, 3)rNsl,N^AX^ 

1,3607179 0,P62«502 

fALFAC^'f 2f 3),Nsl,NMAy^ 

0,0000000 0,0000000 

CAl,FAf^/' 3# 3) fNsU^-^Ay) 

.3,20601^4 .1,8379765 


fi7NT^PTS,= 272 
0>963909 

o ', ooooono 

.0 ,950236b 
fllNT^PTS.s 25b 
-0,6103632 
0.0000000 
0',5291620 


TNT, POlf^^T TOTAL 


1676 



DETF>'hT T ON Of- AEHODYNArilC inFluEmCE 


TjTlk s WfCTAmGIJLAW A» s 2 

IDl = ^ 

ID?. a 

wsro s ? 

pp 5 3 

ft^TVPE s 

MM S T 

MftTE s 

SYTY»F !! r> 

N M A X S ^ 

*^onES 5 ^ 

M UKiSv ^-* s r 

L M T M z 1 

L M A y = Sit 

. TJ^AX = 1 P 1 

J J - ^ 

dE'tao s a.ooooo 

PPS s 0,^0000 

MACH s D,BOOOO 

IK!f-'STO s 1 

f 

0,000000 
0,500000 
1 ,^00000 
I !s:DEY E.TA 

1 i 0 , 7 0 7 U 

2 2 - 0,00000 

3 3 - 0 , 7 7 1 1 


a T p T X 


11 - 13-73 



STAWTING 2000 LOOP 


F.TA= 0,7071068 CHlaO, 0000000 
(AtFA(N, I, 1),N«1 #NmAX) 

0‘,85aa37? .0, i8a3266 

(Al FA(N, n f 

r.H3t>^27 .. ,«0*.l75'>l2a 

(ALFA(N, 3| nrNsl.l^MAx) 

I'.0l7^b37 •,0.0916097 

FTAa 0.^7071068 CH I SQ , 4999998 
C ALFA 1 » 2) f Nsi f ^"'AXl 

Q',97ai887 0,2500034. 

(ALFfcf^, 2, 

2‘, Wl979<i6 0,7071 1 13 

(ALFA.C^» 3^ 2).Nsl,wpAy) 

2',48fl8277 0,7500026 

FTAs 0,7071068 C P I * 1 , 0 0 0 0 0 0 0 

f ALFA (N, 1 , 31 » f^ = 1 • ) 

l', 1422063 0,6843331 

( ALFA (N, 2f 3) » Nsl , N^Ay ) 

2’,9322974 I ,5901354 

(ALFA(N, 3# 3)^Ns1>M^Ax) 

3‘,0470359 1,5916157 

FTAb-O, 0000000 CHlsO, 0000000 
( ALFA C^* I ^ 1 ) • ^sl ' N?LA,X^ 

0,7352671 .0,2561758 

(ALFA(N, 2i 

.0,0000000 0., OQOOOOO 

fALFA(N, 3, UiNsUN^-'Ay) 

.1,1952249 0,1642214 


«TNT,PT8,s 344 
.0,2398086 
•0,2858456 
.0,1686799 
«T^T,PTS,5 408 
0'.3921047 
O’ 7525932 

0'.70Q3559 
ttlNT,PTS,s 344 
-0,2398084 
.0,2858455 
. 0,1686796 
tt!MT,PTS.z 336 
.0,3469a62 
0,0000000 
0‘.27759ft5 



F'fAs.O.OOOOOOO 

) 9 ?) ?Nxl,N^‘A)r) 

ll02BaS6b 0,2500033 

f A I F A C ^ d.) 

o’, 0000000 0,0000000 

f AUF A f *^1 5* 2) f 1 # 5 

.2,5a33286 ,0,7500028 


s?‘‘jT,PTS,s aoO 
0>638i^S2 
0 0 0 0 0 0 0 0 
- 0,7751673 


FT Ar«0 , 0000000 C h I s 1 , 0 0 0 0 0 0 0 iaTNT,PTS,s 320 
(ALFA(M, 1, 

r,l6760^9 0, 7^61666 -0,3^69412 

fAL.FA(W» 2» 31 

o’, OOOOOOO 0,0000000 ■ O.OOOftOOO 

(ALFACN. it 3) 

-3,0926596 -1,6642203 0‘,?775957 


TN'T. PD!WT total 


2l5? 



0ET£R15:J^aTI 0N- , OE...AERaD¥NAMlC... 


title s rectangular '.^ING AR a 2 11 - 13-73 


IDI 

X 

a . 

TD? 

X 

2 

MSEC 

a 

1 

PP 

Ml 

. n 

5 . - ■ - - 

CWTYPE 

X 

0 

HP 

X 

1 1 

mREE 

. X 

n 

SwTVPE 

X 

0 

NPAx 

X 


KK..„. ... 

X. 

-n 

MODES 

X 

1 

UNSyv 

X 

0 

LRIN 

: -« 

.1 

1 . MAX 

X 

511 

jjmax 

X 

191 


.s 



OELTAO 

m 

4,00000 

EPS 

B 

0,50000 

MACh 


0,00000 

infsto 

B 

0 


(CHICP) 

- 0,00 0 . 0-0 0 
0 , ia 6 iia 7 
0,500000 

0. 85.35 5 3- 

1,000000 

I Index eta 

1 I 0,9659-3 

2 2 0,86603 

3 3 0,7071 1 

a 4 , 0,50000 

5 5 0,25882 

6 6 -. 0,00000 

7 7 ^0,25882 

8 8 -0,50000 

9 9 -0,70711 

10 JO -0,86603 

U 11 -0,96593 



STARTING 20'0n m LOOP 


RG^'S 

1 » 5 

have 

BEEN 

STORED 

RO^S 

fe-1 0 

have 

B fc. E N 

STORED 

ROWS 

11-lS 

Have 

BEEri 

STORED 

RO^^S 

o 

!\) 

1 

have 

«tEN 

STORED 

RO'^S 


Have 

BEEN 

STORED 

RC'^S 

2b-30 

H A V E 

been 

STORED 

TNT. 

POT NT 

TOTaI 

5 

BRR 1 ? 



APPEMDIX-I I 


PROGRAM LIST! MG 


55 



li (U u 


DQUBie ETiio PIpSTHETAp AUP A 

PP^e.‘=^'fVPEpgW?Y^EpTnLEeUt^SV(^pU5()y&cU?pUaOpUnpUJia 
XNTESER ODISf^p UCO 
LOGICAL CONTlNpSUPRSSpNOeel pM|fn 

logical FReSHpCONVRSpBaTCM 

real IRREGE, IRREGipH&CH 

RElU IK 

REaL JJSp MODSp MMgp HRESp NMaS^ L^ASp LMISp WSTSp 

U^eSp TWSSp HftCSp NOHSp 
OIMEMSION TITLEfabU CHIcPf205p N|NDEKf^T3 

OlHENslON ET&n®^5o STHETA(58iS3 P HSILIPf38^)p CORDXpeSGS) 
dIhEnSION TLEl(563 3<iTLFR(3a3)pTTEL(S03 3pTTER(SB5 3 
,^^,,0IHEN5I0N IRREGE at least as great as 

.••.fOiMEMSiON iRRfGi at Least as gr^at as nma«*pPo 

• f ^dimension ALFA AT LEAST AS GREAT AS NMA»*KK*PPg 

,,,,,OIMEMSION ALFA? AT LEAST AS GREAT AS NMaX*KK*PP 

D IRENES ION IRREgE( UOO) p IR»EgI C2P5) 

DIMENSION ALFA(aOOO) 

DIMENSION ALFa2(2000) 
dimension IKfa7) 

DIMENSION SKfA?) 

dimension GNL(IO) 
dimension GPNLf In) 
dimension GNT(IO) 
dimension dHNDX f 200) 

DIMENSION HNV2R0f200) 

dimension siNrSfuS) 

dimension fSINE$(A8) 

dimension l..TOTRf205 

dimension NINTP(?0^ 

dimension G(5l20)p XPf512)p HHN(lO) 

dimension ETACPf«7) 

COMMON PI 

COMMON LMAXpUPwRpLMINpUMAXl ,G, XP 

COMMON KpYSpOELTApHHNpNMAX.SUPRESpNOCCIpRLLMIN^USpUTOTAL 
equivalence (ALFA2p IRReGE^ETaCP} 

DATA USp U7pU20pUU/5p 7, 9,11/ . UCO /h/ 

DATA COM /tiU / 



o o r» 


D*TA WINS / 5HWIN / 

D*T* CONS / SHCON / 

D*T* STOS / JhSTO / 

DATS DELS / 3HDEL t 

DITS JJ S / JHJJ / 

0*T* MOOS / JHHOO / 

OAT* PP S / SHPP t 

DATA CWTS / JMCWT / 

DATA mm S / SHMM / 

DATA MHIS / ShMRE / 

data SWTs / JMSWT / 

DATA NMAS / sNNTt* / 

DATA CMAS / ShUMa / 

DATA tMIS / 3HUMI / 

data KK $ / 5HKK / 

Data nsts / jhnst / 

data imps / SHIWr / 

DATA IWCS / JHJWC / 

„ data iws* / SMIWS / 

« Data macs / JhMac / 

data sups / JHSUP / 

data cell / 3 HCC 1 / 

Data nohs f 3hnoh / 

DATA NINS / 3hnIn / 

DATA STAS / 3HsTa / 

OATa TSSS t 3MTSS t 

PI«3', »ai5T2S535e97<»3 On 
U6 a UCO 

C . , ,CONVR8,.TRUE, FOR CONVERSATIONAL USE', 
c!!!,. False, for batch'. 

wRjTr (uSf 20Sa) 
rEaD(USiS) Batch 
COnVRS, not, BATCH 
an CONTINUE 

IF(rONVRS) WRITEfUcO»2o2e) 

READtUSiS) lOl 

aIMFIL is for aHES' tss version onlv. 


SET CONVRS TO 


aImfil issues doef commands to tss. 



call GtMFlLdOn 

REWIND U7 

READfUT) rOt #PP#MM, cWTYPE,SWTYPE,UNSYM,NOL,NDT,MREF, JJMaX,I,TITLE,N^ 

ITITL 

READfU?) (CH1CP( I ) I I«1 rPP) I (NI nDEX (I) #I« liMM) # 

J (TLtlCJ)# J»l» JJMAX)# (TLeR(J)# 

1(TTEL(J)fJ* 1^ JJMAX), (TTEM(J),Jal,JjMAX^ I 

I (ETA(J), (8THETA(J),J*lrJJMAX)f • 

1 (XSILIP(J)* J«t#JJMAX), (C0PDIP(J),J»1.JJMAX), 

1DUM,DUM, aR,OUM,PMaCH 
ENDFILE U7 
JS CONTINUE 
NmAXsPP 
KK*MM 
EPS*, 5 

OElTA0pA«**2 

JJaJjMAX 

MODES«l 

N S E 0 9 0 

LMAXsO 

LMIN.o 

t^PSTOoO 

IWRiTEcO 

lWCCP«e 

IWStPift 

CONT 1 N 9 CONVRS 

SUPRESe, false, 

NOCCi«,TRUE, 

NONE AO* (5 
N1NT«, false, 

JJMAXlsJJHAX^i 

JJUJJHAXI 

MREF1*MREF*1 

MACH»RMACH 

MCALCai 

5‘ft continue 

IE(CONVRS) WRlTEfUcO^S) 
reaO{USiU com 

IF (COM ,EQ, wins) go to 20 



IF 

(COM 

,EO, 

CON$) 

GO 

TO 

00 

IF 

(COM 


STO$) 

GO 

TO 

100 

IF 

(COM 

,EQ. 

DELS) 

GO 

TO 

no 

IF 

(COm 

• eo. 

EPSi) 

GO 

To 

ns 

IF 

(COM 

• eo, 

JJ $) 

GO 

TO 

120 

IF 

(COM 

.to. 

mods) 

60 

TO 

130 

IF 

(COM 


PP %) 

GO 

TO 

lao 

IP 

(COM 

,fq. 

C*^TS) 

60 

TO 

ISO 

IF 

(COM 

.EQf 

MM $) 

60 

TO 

160 

IF 

(COM 

.EO, 

MRES) 

GO 

TO 

ITO 

IF 

(COM 

• EQ, 

SWT$) 

60 

TO 

100 

IF 

(COM 

i^O# 

NHAS) 

GO 

TO 

190 

IF 

(COM 

.EO. 

LMA.S) 

60 

TO 

200 

IF 

(COM 

.to. 

LMIS) 

60 

TO 

210 

IF 

(COM 

,EQ. 

KK S) 

GO 

TO 

220 

IF 

(COM 

,EO, 

NSTs) 

GO 

TO 

230 

IF 

(COM 

,EO, 

IWRS) 

GO 

TO 

2^0 

IF 

(COM 

.EQ. 

I^CS) 

GO 

TO 

2S0 

IF 

(COM 

.EO, 

IWSS) 

Go 

TO 

260 

IF 

(COM 

.EQ. 

MACS) 

GO 

TO 

270 

IF 

(COM 

.eq. 

SUPS) 

GO 

TO 

200 

IF 

(COM 

.En. 

cds) 

GO 

TO 

290 

IF 

(COM 

.EO. 

NOHS) 

60 

TO 

295 

IF 

(COM 

.EO. 

NINS) 

GO 

TO 

296 

IF 

(COM 

.£0. 

STAS) 

60 

TO 

300 


F0« aMESi TSSS version ONLY, 


IF (COM ,EQ, TSS|) CO TO fcO 
IF(j^ONvRS)NRlTE(U5,ao23)cOM 

IF(CONVRS) GO TO SO 
WRITE(U6,2023) COM 
STOP 13 
CONTINUE 

REAO(Us? 7) (ETACP(n#I»l# 20 ) 
CAUL OBEY (00# eTACP) 

GO TO So 

contin«,not,contin 
GO TO SO 


ALOWS T58 command TO BE ENTERED 


too CONTINUE 
STOP 

lift READ )DEtT*0 

GO TO 5ft 

Il5 READ(U5r« )EpS 
GO TO *50 

120 REA0fUS,5 )JJ 

JJsMiN0( JU,UJMAX) 

UJ 1 sJiJi 1 

00 122 UI.JJMAx 

iFCNODf JJM aXi , JJ i ) ,t0,0) GO TO 

122 JJlsJJUl 

123 JJaJjUl 
GO TO 5ft 

I3ft REA0 cU5,5 )N0DES 
GO to 50 

UO REaO(U5?5 )PP 
NmAX«PP 
GO To 5ft 

150 REaD(U5,5)CWTYpE 
lFfCNTYPEM5li I52fl57 

151 READ(U5#a)(CHlCPfl)»I#l,PR5 
GO To 5ft 

152 lFfPP*2) 156, 155, 153 

153 IMaXcPP,! 

DTHETAaRl/FLOATtPP-i > 

DO 15^ I*2,IMAX 

l$tt CHlC:Ptn«(l.-C08(rLQATtI.n*DTHETAn/2 

155 CMICP(PP)«l, 

156 CHICP(l)«0 
GO TO 5ft 

157 DTHETA«PI/FLOATTPP+PP+n 
00 158 1*1, PP 

158 CH!CP(n»n.*C0SfFL0ATcl + n*0THETAn/2 
GO TO 5ft 

I6ft READfU5,5 )MM 
KK«MM 
■MRpFaMM 
M«EFl»HReF+l 



GO to 50 

ITO REAOcUSfS) href 

M»EF«HAXO(lfMREF) 

HREF«HINO{MREF, JJ) 

MREFiaHREF+l 
DO 175 

IF(HODfsMl»HREFn,r0.0l GO TO 175 
173 HREFl«MREF!.l 
175 MREF»HREFUl 
GO TO SO 

l$0 READ(U5#5 ISRTYPE 

IFCSWTYPE) l«2rl«5fl82 
182 REAOCUSfS >(NiNOEXcn»Itl»HM) 

GO TO 50 

1«5 DO 187^I«i,MM 
187 NINoEXfI)«l 
GO TO 50 

IQO «£A0(U5i5 5NMAX 
60 TO 50 

aOO READ(U5rS )IHAX 
GO TO 50 

2lO READfU5f5 )IM1N 
GO TO 50 

220 REA0(U5f5 )KK 
GO TO 50 

230 INFSTO»HOO(lNFSTO^i f 2) 

60 TO 50 

2ilO IWRTTE.MnDf lWRlTFflj2l 
GO TO 50 

250 IWCCP«HOOf 1HCCP^||2) 

60 to So 

280 n«SCP#MoO(lH5CP+lf2) 

60 TO 50 

270 RfAD(U5#il )HaCH 
GO TO 50 

280 SUPRFSs,fgOT,SUP«FS 
GO TO 50 

2R0 NOCCI8.NOT.HOCCI 
GO TO 50 



295 NOHEAOaMOOfNOHEAO^ W2) 
r*0 TO 50 

296 NINT«,N0T,N1NT 
r.O TO 50 

500 CONTINUE 
nSEQcNSEQ^} 

ID2«0 

IFflNFSTO.Nr ,0)GO TO 505 
C 

C,,,,,FOH aNE8( T88 SYSTEN ONLY, DEFINES OSNaME OF IDFILE. 

CALL OBEY(22,2<iHDOEF F T09F00 1 m 1 OF ILE ) 

C 

rewind U20 

REaO(U 20) 10| 102«ID5»IDa 
I02»l024t 

rewind 020 

WRITE (020) lDrl02»lD5»lDa 
ENoFILE u«o 

^ Call osEynb, i6hrEi ease fto9fooi > 

<=» IF (NOHEaD.NE.O) WRITE(UCO»2027) ID2 

505 CONTINUE 

iF(CONVRS) WRITE(OcO,2o29) 

«EA0(US>5) 00I5W 
IF(OOlSH,EQ,0) GO TO 50a 
U6su 

OOISK«MOO(MAYon»ODlSK5 1 JO) 

W«ITf(Ocnf2026) ODISK 

CALL obey f 16, ISHRELEaSE FTO^iFOO) 1 

CALL CVRT(00I8K,1, 

1 aaHCtODEF FToaFOOW,OOTRUT,UM,Nt, jl,6JO 
2ETACPi8,6H(SAa) ) 

CALL OBEY (52fETACF) 

rewind 06 

IF(cONVRS)WRITE(UCO,2o2T) ID2 

5o« continue 

Call inital(Cwtype,pp) 

IF(NOHEaO,NE,0) go to ^10 
WRITE(U6,2000) 

WRITE (06,2001) (TITLE(N),Ntl,NTlTL) 



rt f>rt KM noon 


wRITE(U 6.2002) I01.lD2fNSEQ 
wRITEfU6,i00S)PRiCwTYPE 
WRlTEtU6,ZftO«) MM, MR£F, SWTVPj 
MRITE(U 6,2005) NM*X, MK, MfiOES, UNSYM 
WfilTE(Ufc,2006) LMIN, UMAX 
MRITe<Ub,?007) JJMaX,JJ 
WRlTE(U6,2008> OeuTIOfEPSjMACH 
WRITE(U6,200<») iNfSTO 

IFflMCCP.NE.e) MRITE(U6,2020) f CHlCP ( I) , I » 1 , PP ) 

iPawSrp.Eo.O) r.o TO Jio 

JRATIO.JJMAXI/MREPI 

WRITp(U6,2025) 

DO 505 111, mm 

J«NINOEX(n*JRATIO 

ETASiETAtJ) 

50S MR1Te(U 6,2022) I,NIMr)EX(I>,ETAS 
310 CONTINUE 

If (MOO(JJM aXJ, JJ l) ,NE,0 )PaUSE'2' 

If (MOO(JJI|M«|fn ,N e, 8 ) PAUSf ' 5< 

IFflNfSTO ,Nf,0) CO TO 1«0 

computing SPaNWISE control points for aIM file; the 

.•.•.ONLY place ETaCP is used is in writing the introductory 
,, ...RECORD of The AIM FILE’. 

JRATIO»JJMAXl/MRFf 1 
DO Sj2 lei , mm 
J.NINO fXf DAJRATIO 
ETACP(i)«ETA(J) 

12 CONTINUE 

FOR aMESi TSS version ONLY, aIMFIL ISSUES OOEF COMMANDS TO 

THE TSS operating SYSTEM'. 

Call AIMFILf I0l»l02) 

NRITEfUl I ) lOl ,1D2,NSED, TITLE, NTITU,PP,CRTYPe,MM,MREF,SwTYPE,NMAX 
t KK, MOOES, UNSYM,LMlN, UMAX, JJMAX,JJ,M aCh, EPS, CELT aO, 

1(ChICP(1)«I«1»PP)» (NiNDEx(I) • isl »MM) , 

1 (ETACPf I) , lil ,MMJ 
3tt0 CONTINUE 



Call INFMATtHM^MRpFrf'P^NMAXWjUNSYH^KK^HO^tS^JJrMCALC^JJMAX^INFSTO,, 
II'*«ITE,U5,UU ,HACH,DELUO,ETA,STHETA,XSlLlP,CO«OrJ^,CHlCP,^^INDEX,Nl. 
INT, lR»Ef,E, IRREG3 ,ALFA,ALf*2# 1K#S«#GNL,(;»»NL,(5NT,OHNdX,hNYZRO, 
l«Hl, J[jMAXl,STNES,CSiNtS,LTOTP,NiNTP,ERS,TLELf TLER^TTELfTTfcR) 

E^^nFILE Ltli 
GO TO 5o 
j FOHMAtfAS) 

% rOkMATf^H ♦ ) 

a FO«HAT(er to.o) 

5 FORMATf^felS ) 

6 FORMAtfLl ) 

7 FORMAT(20Aa) 

2000 FOKMiTf JHi/a7N0DtTE«HlNATI0N OF AERODYNAMIC INFLUENCE HaTWIX //) 

2001 FORMATf^HOTlTLE^"X>^N« #?0Aa^ 

2002 format fflHOIDl #6X, IH«, IS/aH I D2 , 6X , 1 H» , I5/SH n 5F 0 f b X ^ 1 H* , 1 5 ^ 

2003 FORMaT(3H PP,7X» |H»,IS/7H C «T YpE , JX * j h« , I S ^ 

200a FOkMATf3H MH,7Xf iHSf I5/5 h M»EF , 5X , 1 h* , ! 5/7h SwTYPE ^ 3X , 1 h« » ^ 

2005 FORMATfSH NMax i SX , 1 Ho , 1 5 /3H KK, 7X, IH3, 15/6H MOOES , a X , 1 H« , 1 UN8. 

lYM^ajr, I51 

2006 FOPMATfSH LMl N , by , 1 Ha , 1 5/5M LMAx#5x, 1H», IS^ 

2007 FORMATfOH J JMAX , ax i I H« , IS / 3 h J J , 7 X f 1 M# , I b ) 

2006 F0RMAT(7H Df LTA0,5 X, IHs^FM , b/an fPS, 6X, ,b/ 

IbH MACH^SX» 1 Hb,fU ,b^ 

200R F0RMATf7H I NFST 0 , 3X M H* ^ I b 1 

2020 FORMATfBH t CH IC P ) / H X # F j 0 , 6) ) 

2021 FORMATMRH j jNDEX ETA ) 

2022 FORMATflX,14,I5#F10.S) 

2023 F0RMATf2aH UNRECOGNISED COMMaNOi 
202a FORmaTmSH EnTER BaTCh ) 

2025 F0RMATM5H ENTER OOlSK ) 

2026 F0RHAT(2RH OUTPUT IS ON , , , OUTPUT . a I M. N ,H,3h,.,^ 

2027 ^ORMATfbH I02«, l3) 

202a FORHaTuSN enter 101 ) 

END 

aUpROUTTNt INPMaT(HM,«REF^PP,NMaXR,UNSVM,KK,H00E8,JJ,MCALC,*IJMAX, I 

INF5T0.IRRITE,U5,U1 1,MaCH,DELTaO,ETA^STH£Ta,XS1L1P,COROIP,CH1CR,NIN 
1DEx#NINT,I«REG€,I»«EGI, ALFA,ALFa2,IK, SX, GNL, CPNL, GNT, 



IOHNDX, HNYIRO^ JJMaXI, sines, CS!NES,LT0TP, NXNTP,£PS, 

lTlEl.f TLE«rTTEL,TTER) 

double PWECISION ETa,PI,SThET&,S1wES,C8InES 
DOUBLE PPECISION DK,0K2,DETA,D£TA2,8JM|P4Jl,y^AirA,RJJja 
INTEGER PP,c«TyPE,8NTYPE,UN5YM,U5pU6,P#PL,PU,TjRAT0,TJJMXl,TjXlMl, 
lun,. GTOTAl 

LOGICAL $YH,ASYH,TLE,TTE,TlNT,N0Tl,N0TT,N0Tl,MNzrR0,NEl,8UPRES,N0c 
IC1,NINT 

real Ik, IR»ECE,IRREGI,MaCH, jOLM,jOLPf JUM, JlLP, JITM,j1TP 
DIMENSION NINdEK(«m) 

DIMENSION CH1CP<RR) 
dimension eta^jjmaki) 

DIMENSION STHETa r JJM aKI ) 
dimension XSiLlPf JJM aX) 

DIMENSION cOROlP( JJM aX) 

DIMENSION tLEL(JJHa), TLERfJJMAX), TTELfJJMAX) 

DIMENSION TTE«(JJMaX1 
dimension IRREGEf2#3>NMAxR#KK) 
dimension lRREGlfNMAXR,PP) 
dimension IK(MM) 
dimension SKfMM) 
dimension aLFA(NMAXR,KK,PP^ 

dimension ALFA?(NMAXR,KK,PP) 

DIMENSION GNL(NMaXR) 
dimension GPNL(NMaxR) 

DIMENSION QNTfNMiXR) 

DIMENSION OHNDX(NMaXR#PP) 
dimension HNVlROfNMAXR»PP) 
dimension 8INES(MM1) 

DIMENSION cSINESfMM) ) 

dimension lTOTP(PPI 

dimension NJNTPfPP) 

dimension G(5120), XP(S12), HHN(IO) 

common pi 

COMMON LMaX,LPwR,LMIN,LMaXI,G,XP 

common X,YS,OELTA,HHN,NMAX,8UPRES,NOCCI,RLLMIN,Uf iLTOTaL 
data CI/>5661R772S/ 

DATA C2/. ,1185098061/ 

DATA 05/1,275239505/ 



O O iJ 


DATA Ca/*,6J66197T23/ 
data C5/1, 3353335/ 

DATA C6/,a50lS81579/ 

DATA CT/., 30010^4386/ 
data C8/J25/ 

DATA C9/,l ,380294361/ 

DATA CIO/-, 8333333 Ul/ 

DATA Cll/-,0625/ 

DATA Cl2/-,25/ 

MRATlO.JJMAJd/fMRjTF^n 
JRATlOiJJMAXl/JJl 
rjRATOaJRATIO^JRATlO 
FTA(JJHAXn««l ,00 
STHETAf JJMAX1)«0,D0 

C..*. *070X41 NlLl BE ThI grand TOTaI OF THE NUHi^ER OF 
C,*,.*1 NTEg«aX10N POINTS USED FOR aU TmE CONTROL POyNTS 
GTOTAL »0 

IF (UNSYMJ llO.lOOiUO 
100 SVHB.TRUp, 

MMPt(HM*n/2 

IF (MOOES, LT,0) MMPbMm/2 
GO TO 120 
110 SVH *,FAlSE, 

MMPtMM 

120 ASYMf ,NOT,SYM 

KU2«MAX0(3»KK*1 ) 

TJJHXUJJMAXI^JJMAXI 
TJXlHl.TJJMXl.l 

.*♦. MOOES, lt,o***odo modes 

,,,,mOOE8.EO,0,,,EvEN wOoES 
,,,.mOOE$,GT,0.,.aLL MOOES 
lF(MOr>ESH30il40,l*50 
130 Kt*2 

KJUMPb^ 

GO TO 160 
140 Kl«l 

KJUMPe? 

GO TO 160 



ISO Kill 

KJUHPci 
J60 rONTJNUK 

MNZERO«MACH,NE,0, 

BETAti; 

ir fMNZERo) BET A ssort ( 1 ,-MACH*»ai 

TSETAsS.sBETA 

BETA?«BETA*«? 

NflsNRAX.rO, 1 
TL£«CHICP(n.LT,l,E-fc 
TTEsCHlCPfPP) ,GT, ,>999999 

TINTb(PP,gT, 2) ,OR, (PP.EU.2 , ANp, ( ,NOT , TEE , OR , .NOT , TTg ) ) , OR, .NOT, 
1 (TLE.OR.TTE) 

NOTLs.NOT.TLE 

NOTIs, NOT, TINT 

nOTTs.nOT.TTF 

PL«1 

PUsPP 

«. lFi;TLt)P(.s? 

>^ 1F(TTe)PUsPP-1 

COMPUTING quantities ASSOCIATED wITM CONTROL POINTS 

C ON leading anO trailing EDGES (EOS', 52,SJ,S6,AND sT) 

Sf,Nsl , 

GNLtllsCi 

GPNLtniCE 

GNTfJjsCl 

if(neI) go To aso 
Do 200 ns2,nhAa 
GNL (N)bO 

GPNL(N)bCJ*FLOAT(N.1) 

GNT(N)sS1GN(GPNL(N),SGN) 

SqNsvSgN 
200 CONTINUE 
250 CONTINUE 

IF(N0TI)G0 TO SoS 

Ci... .COMPUTING quantities ASSOCIATED WITH INTERIOR CONTROL 

POINTS (EOS, 6S AND 6«5 

00 JOO PbPl,PU 
*s2,*CmICP(P)-1 . 



n.o 


, pmi«arc6s(«x) 

SPHi*5QRT ( I 

OHNOX( 1 f P)»Ca*n ,*X) /SPHI**5 
lF(NrnGO TO 500 
DO 290 W«2 ,NmAx 
0UMl®N«l 

DHNDX(NfPj*Cl*0UMl*C08fDUMl#PHn/8PHl 
290 CONTINUl 
300 CONTINUE 
305 CONTINUE 

RJJl<*i,25D0*«JsIl 

OHcON«OClTAO/PJJJtii*2 

C,.,, •DETERMINING MN OF frQ\ 1 FOR THE CASE Y«0'. 

00 550 P»1,PP 
X»2,*CHlCPfP)»l, 

Call infzY(X#nmaxihhn) 

00 55 O Ntl jNmAx 
4T 350 ■hny1RO<n#R)«hhn(n) 

^ WRITE (U8^ n 

c****************** 

bO 2000 M«McALC#NHP 
M2»NINPEX(H)#MRATI0 

.♦.fUTOTPeP ). number of integration points used over the 
.•••ENTjRE' SPAN FOR Each ChOROwiSE CONTROL POjNT location, 

DO 560 P«I,PP 
LTOTP(P)tO 
360 CONTINUE 

IF ifNlNT) W«ITE(U(!!ii7)M 
BB2C«BETA/C0RDIPfM?) 

Bat«6B2C;B62C 

FTAMS»EtA(MJ) 

IF (NOTL) GO TO 500 

c computing quantities associated with leading edge 

Cf*«.#CONTROL POINTS AND WHICH VARY AS THE SPaNWISE POSITION VARIES, 
TANLEL«TLEL(M2) 

C08PHI «B|T A/SORT (BET A24TaNL£L**?1 
EMs,5*( I ,ATANLEL*COSPmI ) 

emcomp*i,«en 



OUMlaSOftT(662Ci^eoSPHI) 

DUMtS*DUNl**S 

C,....COS, 7<t AND 7S 

jouH,«,*ouM|*cELJ(e«,eM,o,) 

JtLMac5*DUNtS*CtU2(EH,EMr2,*EM*EMc°PP) 

TANLE»aTLE«(HE) 

tP(AB8<tA*ILELAtANLERt,LT.l,E.6l 60 70 <100 
COSPHIaBrTA/SOPT (B eTA2^TanCeR**E) 

EPa,5*tl,*TANLC«*e08PHl) 

EHeOMPat,»EH 

DUMiaSaRT(BB2C/C0SPHl) 

OUH|5aOUMi*a5 
GO to «10 
«0n tOMTtNUE 
DOH?aEM 
EMaEHCOMP 
EMCOMPaOUNj» 

«10 CONTINUE 

J6XrB«,ADUHiaclL2»e«»tP»0 ) 
JJLPaC«*0UMJJaCEl,2fEP»EP*2.*EP*E*<C0HPi 
500 CONTINUE 

IFfNOTl.AW.NDTTi 60 TO TOO 

EO. 00 

CMla(l,«Ef AMs)*,5 

C,.,.,EMIcON ISTHI eOHPLlHENTARr PaRaHETE# TO EmT 

EMtcOHa(l.»ETAM8)a,9 

DUHtai,*EMlaEMtCOH 

C EOS, SB.5B. ALSO SEE EOS, 70.77, 

GAMJNtC0*CEL2(EMl, I , , l,.2,aFMl ) 

60MtMse7*CEL2(ENl«ENi,()UNt7 
e*HSPiefr*CEL2 (¥N I eON , I , , 1 , -2 . aEn I cOn ) 
6AMtPaC7*CEL2(EMtCOH,EM|cON,OUM| ) 

C.,.,,C0NP0T1M6 OUaNTITIES aSSOCIaTEO WITH TRaIUnC EOSE control 
C..,,, POINTS and HhICH VaRT a 8 TmE SPaNNxSE POSiTiON VaRTES. 
TONTELaTTELtM*) 

C0SPHlaBCTA/S0NT(BfTA2«TANTEL*«>) 
iNa,5A{l',sr4NTELice»NNr) ' 

EMCOMPat,.EH 

DUNl aSORT (BB2C/C0SPH1 ) aa J 


68 

JlTH»C5<»OUMi*CEL2(EM,EM#a,*EM*EHCOMP) 

TANTERsTTE«(m?) 

IF (4BS{TANTEl+TANTtRi,LT, 1.E.61 GO TO 600 
CCSPHUBFTA/SQRTfBf TA2 *TaNTf.R** 2^ 
f I ,«T ANTfeRttCOSPHl ^ 

OUMltSQHT (B82C/C08PHI)**5 
00 TO fein 
600 CONTINUE 
OUM?«EH 
EM.EMCOMP 

Emcomp»oum? 

6ln CONTINUE 

Jl TP#C5«^OUMi*cEL2fEM,EM,2^AEM*EMrOMP^ 

700 CONTINUE 

C,,,,,C0HPUTING quantities associated NITh interior rONTROL POINTS. 

SINES and cosines IN EG, 22*2S 

^ J«r(KL-n*H2^TjXlMl 

00 00 flOO K#KL,KU2 

INdEx«HOO( JJMA y 1 ) ♦! 

If (MOOf Ji TJUMXD + l-INOf XlTlO# 7l0i720 
710 5INES(K)*8THFTA( index) 

CSINES(K)bETA(INdEX) 
r,o TO 75n 

720 SlNES(K)e-sTHETAclNpEX) 

CSlNES(K)tt«tTA(lNOEX) 

750 J«J+N2 

c. . . t *En, 10 

eoo SK CK)oSTHETA(H2)tCSINEScK) 

C computing IK OF fOS, 22. 2S 

DO 900 K#KL#KK,KJUMP 
IF fK-2) 8l0^820f 650 
610 tKn)fcC6*(CSlNEsf5)^C9) 

GO TO QOO 

820 1K(2) *n0*r2.D0*SlNES(5)*SlNES(2)TcSINES(5)*CSINES(?)) 

GO TO 900 

830 IF fK«.5)eao,e<io#850 

6<I0 lK(3)*CU*f USDO*SINFS(4)*SjNES(2) + ,SOO<»CSINESra)»CSlNESf2)fCSINES 


1 (3)^C9) 

GO TO 900 
850 DK»K 

DKinKmB 

IK(K)sCl2*({OK*8INt'S(K + n*$INE:S(?)*CSINES(KfUtc8INES(2))/COK*OK-i 
1,00) • (0K2*8INES(K«n*SINE8(2)*r8INES(K.l)#rSINES(2) )/(DK2*0K2- 

U.DO)) 

900 CONTlNUe 
87 

BJH«RjJia/STHETACM2) 

tt,#fCOWPUTING the exact integrals OF THg IRREGULAR 

.Parts and the contribution to alfa fpom etai.eta, 

IF (NOTL) r,0 to 1000 
OUHT80'. 

DUM8#0', 

IF (TaNLEL.NE, T anted 0UM7».^5ik(TANLEL«T anted /CORD lPfH^) 

If (TaNlER,NE,TANTER) 0UM8«-,5fk(TANLER«TANTER)/C0R0lPfH2) 

00 950 KiKL,KK,KjUMP 
OUMl«JoLM*CSlNESfK) 

DUM2#J0LP*c8INESfK) 

DUH6aFU0AT(K«n#SlNES(K)/STHETA(M2) 

DUH3«JoIH*DUH 6 
0UHR#JiLMAf8INBS(H) , 

DUM5«40LP*DUM6 

DUM6«JiLP*CSINF8(K) 

OO 950 N«I,NmAX 

IRREGEf W lf»^#X)*GNL(N)*OUMi 

IRREGE(2M»N,K5«GNL(N)*0UH? 

IRREGEf l»2^N,K)»*GNLrN)*DUM3tGFNLfN)ADUMaTDUM7<rlHREGECl»l>^»K) 
IRREGE(2,2 ,N,K)«GNL (N) aOUHS^GPNL (N) aDUM6 + 0UM8*IHREGE r2, UN,K) 

EG. 51 

AUFAfN^KD)«IPPEGE(l,l,N,K)^GAH3M*IRREG£f2,l,N,K)*GAH5P*IRREGEn,2 
I »N#K)*GAMlMflRREGE(2,2 » N,k)*GaMIPt8jM*8k fK ) ♦HNYZH0(N» 1 ) 

950 CONTINUE 

1000 IF<N0TI) GO TO 1100 
BBC2 • BBC*a 2 
DO 1050 P*Pl,PU 
00 1050 Nii.NmAY 
IRReGI fN,P)B,5*BBC?*DHNDY{N,P) 


DO tObO KeKl,KK,KJUMP 
CO, 2b 

AlFAfN,K,P),I«RfcGX fN,P)*lKfK)*ftJM*SK(K)*HNYZ«0(N,P) 

1050 CONTINUE 

IlOO iFfNOTT) UO To 1200 

DO U50 K«KL,KK,KJUMP 
DUM^*«CSINESfK)tJ^TH 
OUM2«*CSlNFSfK)*JlTP 
DO 1150 N»1,NMAX 
IRPCGEf 1 ' 

IRReCE(2f 3rN,K)»DUW<>*fiNT fN) 

C,,.,«CO, bS 

aIPA(N,K,PP),ZRReGE( l,3,N,K)*GAHlM4lR»EG£(2p5^N^K)*G*MlP*8J«*SK(X) 
l*HNYZRO(NpPP) 

1150 CONTINUE 
1200 CONTINUE 

JU*MODfM;?, TJRATOI+JRaTIO 

start Of SPaNWISf INTfGRaTIOn **«* 

00 IbOO J»JL,JJMAX,TJRat6 
C..,,, determining SK(ETaOi EQ, t«» 

L»fKL-l)*J+TjXlMj 
LJUMP»J^KJUMP 
00 I3bn K«KU,KK,KJljMP 
INOfx«MOD(Li JJMAY l>il 
SK(K)bETA( INDEX) #STH£TA(J) 

IFfMODd# T JJMXj ,GT, INoEX) SK fK)».SKfK) 

13<>0 L»L*LJUMP 
C*,f,,EO,b0 

0ETA»ETA(M2)«ETA( J) 

OETA2»oETA*oETA 

BJM,-5TNETA(J)/(«JJI*0ETa2) 

Y«TBETA/CORDlPrJ)*OETA 

YS«Y*V 

C LIMITING THE MINIMUM NUMpER OF CHOROi*iISE INTEGRATION 

c stations, 

RLLMlN»,ObeRA/(Y$t; 00755) 

o«*«,*eo, ro 

ABDETA«DABs(DETA) 

OELTA*OEtCON/(SNGLf«BJM))**EPS 


DUMJt2,/C0«DlPf J) 

t>UHl *DUM2* (XSH IP(M2 ) -XSILIP CJ) )«l . 

OUH2#OUM2*COROIPfMa) 

8oETASsS9RTrA0DEtA)*STHETA(J) 

I»1 

IF (DETA,LT,0,00)I»2 
iFfNOTL) CO TO ia?n 
x«dumi 
CALL INFFUN 

LTOTP{t).tTOTAt^LTOTPn) 

IF(NINT) NiNTPf n«lTOTAU 
00 laoo K»KL,KK,KJUMP 
DO 1400 NsI^NmAx 
C 0* QP 62 

DSKHN«(IRREGF(U UN,K)*IRRCGE(I,2#N,K)*A8DETA1*S0ETaS 

Co. 51 

ALFA(N,K. I )«AtFAfN,K, tK)*HHN{N)-DSKHN) 

laOO CONTINUE 
1A20 IF (NOTI) go To 1440 

OUM380FTA2*AL0Cc A eoETA) 

00 I4jft PaPl.PU 
X#DUMUDUN2*CHlCPfP) 
call INFFUN 

tTOTP(P)»LTOTALALTOTPfP) 
lF(NINt)NINTP(P).LTOTAt 
DO 1430 K«KL,KK,KJIUMP 
00 143o N«1 ,NmAX 

14J0 ALF*tN,K#P)*ALF6fN.K^P) ♦ B JM*SK (K )* (HHN f N )♦ IRRpG I f N • P ) *0UH5) 

1440 IF (NOTT) GO TO 1460 

X«DUMj+0UH2 

call INffUN 

UTOTP(PP)»LT0TALFLT0TPfPP) 
iFfNINT) N1NTP(PP)eLT0TAL 
0UH4»ABDETa*SDETaS 
DO 1450 K«KL,KK.KJUMP 
DO 145o Nil.NMAX 
C.....EO. Tl OR T2 

DSKHN.IRREGE ( I » 3 # N , K ) *DLIM4 

1450 ALFAfN,K,PP)»ALFA(N,K|PP) ABjMifSK(K)*HMN(N)«DSKHN) 


ia60 COMTI.NUt 

iPffilNT) WRlTr(U6»8) (nInTP(P) ,Pal ,PP) 

1600 CONTlNUt 

***<k***«^* SPiNNIse integration finished *♦*★*#***♦ 

c 

c ,,,, .changing to MOOES SKgSINfK^THtTA) 

c 

IFfKK.FQ.n GO TO 1609 
Kl^sKl^S 
OO 1608 P81,PP 
00 1608 Nel ,N hAx 

' IF (MOOES.Nl.O) ALFAfN,2,P)«2,*At.FA(N,2,P) 

IF(KK,LT,5) go to 1608 
f)0 1606 K»KL2,KK,KJUHP 

1606 AIF ArN^K,P)a;>,*ALFA(N,K,P) + AUFA(N^K*2,P) 

1608 CONTINUE 

1609 CONTINUE 
C 

C.....ENO OF changing modes 

c 

IF (InFSTO.NF ,0) 60 Tn 1610 

NR0W;»*PP*M 

nR0W1«nR0W2«PP+1 

c 

THE INFLUENCE MaTRIX IS STORED aS SINGLE PRECISION 

C 

00 1605 Pal ,PP 
DO 1605 KaKL,KK,KjljMP 
00 l6o5 Nbi^NmAX 
1605 ALfA2(N,k,P)*ALpA(N^K>P) 

WPlTE(Uin (((ALFA?(N,K,P),Nirl,NMAX),K«KL,KK,KJUMP>,P,l,PP> 
WRITF (U6#5)NR0W^ ,NR0W2 
1610 IF f INRlTfc.^Q.Ol GO TO 1650 
DO 1630 Pal,pp 

WRITE(U6*2) ETaMS,CHICP(P) ,LTOTP(P) 

DO 1650 KbKL,KK,KJUMP 

WRITf fU6»5) K,P 

write (U 6» U) (ALFA(N,K, p) ,N,1 ,NHAX^ 

1630 continue 


r»r»r»r»r»r»n rir»r»or>o 


165P continue 

DO ll<>9 PeJ,PP 

GTOTAU»GTOT*L*1-TOTP{P) 

1160 CONTINUE 
2000 CONTINUE 

WSITE(U6,6)GT0TAt 

RtTURN 

1 format t22HjST»RT INC 2000 N LOOP //) 

2 FORMAT (ShOETA*. FIO.7, 6h Ch1*i P^tT, Um *INT,PTS,« >i5) 
i format (RH (ALPA(N, ,12, IH, , 12, llHwN»I,NMAxn 

U format (iX,Ft<I.T,SFl5,7) 

5 F0RMAT(5h0R0W8 ,1J,Ih*» Hi ITM HAVE BEEN STORED ) 

6 F0RMaT(20M0INT, point total a ,Ifl) 

7 FORMATfRHflM a ,IS) 

B FORMATMX, 26IB) 

End 

subroutine INFFUN 

c THIS subroutine determines the pressure mode influence 
”|*|fUNCTIONS, which 4RE the INTEGRiU ON OF 

,,,,, HN(XP)*XBA«(X.XP,Y)/j, 

,.,,.hN are ThE V»n SPEiCEL ChOROWiSE FUNcTjONS 
..,..THE integration is done REPEaTEOUY with 
.....TWICE AS Many integration points Each time 
..,.,ONTiL either ThE maximum ChaNqE is less 
.....than delta and LL.GE.HLLMIN or until LLeUUMAXl'. 

LMIN is the starting number of integration points, 

.....NESTins is used so ThAT after ThE first iNTEr, R ation ThE 

Values of kbaR need only be determined at 

....•EVERY other station, 

The integral is computed using the quadrature formula 

Of problem S7, P,S6« of HILOEBRaNDiS numerical methods. 

LLMAX a(LMlN*j)*2*»(LPwR-n-i 

double precision sum, Pi 

DIMENSION G(10,5l2),XP(5l2),HHN(10),HHNRf 10),SUM(10 ),dIFFN(10) 
REAL KBAR 

integer U6 

LOGICAL SUPRES.NOCCI 
common Pi 


r» o o ri r> 


COMMON LLNA><,LPW«»UMlN,tLHA)(I YS,OEtTA,HHN,NMA^^SUP»E5,NOC 

lCU^<ULMIN,Uh,LTOTAl 
XMXPeX^l , 

KBANr (1 , + XHXP/SQRT f X mXP*XMXP + YS) )/2, 

DO 10 N*1,NMAX 

SUM (N)aK8AR*G(N, UMAX! ) 

10 HHNfN)tO, 

20 CONTINUE 
Ll»UMTN4t 
LjUMP*ttMAXl/LL 
IJUMP?«UJUMP 

c,,,..each Pass through the j loop represents a chordwise 
c,.,,, integration 

DO 200 Jsl,LPWR 

c,,.. .PREVIOUS results are Saved in hhnp; 

00 So n»i,nmax 

HHNR(N)aHHN(N) 

So CONTINUE 

o,C,.,.,LL» number op integration POINTS INCLUDING T^E 
■^C,.,,,ONE aT THE leading EDGE, 

DO 150 LkLJUMP,LLHAX,LJUMP^ 

XMXPxX»XP(L ) 

KBARtiJ,4*xMxP/$ORTfyMxP*XHXP*YS5 
DO 150 N#1,NhAX 
SUMeN),SUMfN)AKBAR*G('N,L) 

150 continue 
PIL*LL 

LjUMP^cLjUMP 
LJUMP»LJUMP/a 

OlFPxO 

DO 170 N«1,NmAx 

THE following STaTEMENT INSURES THaT ONLY HHN»S COMPUTED 
with a certain decree of precision i^iLL p£ USED IN ThE 
comparison convergence TEST, 

HHNfN) WILL t HHNlNFfN) UNTIL THE REQUIRED D.O.P, IS 
obtained, 

lF(LL«N,LT,«nGO TO l75 
HHN(N)aSUM(N)/RLL 
OlFFN(N)»ABS(HHNR(N)-HHNtN)) 


oooooooooor* 


170 OlFFBiMiXlcOlFfjOlFFNIN)) 

175 LLsLL*LU 

180 CONTINUE ^ _ 

IF (OIFF .IT, delta .and, »LL , GE .HULMIN) OO TO JOO 
200 CONTINUE 

LT0TAU*UL/2 

IF fSUPRE8)RETU«N 
WRITF fU6H) 

GO TO So5 
SCO continue 

LT0TAU»LL/2 " 

IFfNOCClJNETURN 

J05 write (U6,2) X,7S»ftElTA,LTOTAL.{DlFFN(N),N»l,NMA75 

Turn 

1 FORMaTC///' THE CH0R0WI8E INTEGRATION FAILED TO CONVERGE') 

2 eORMATC X »iiFl5,6»i VS t'^FlStO/' dELT* «i»F12,8.' 

IUMBER of integration points *',I5/' DIFF • . , 5F 12 ,8/8* , 5F 1 2, 8 ) 

End 

SUBROUTINE INITaL(CRTYPE.PP) 

THIS subroutine computes ,c,, the Values of the 

....iChORORISE pressure nodes evaluated 

.....*T POSITIONS given by the aRRaY fXP), WHICH IS aLSO COM. 
.....Each mode haS been multiplied by pi/?«8IN{PhI) . 

.....This version is for ThE van spejgel pressure modes'. 

.quantities are determined at EACH XP s -COS(PhI), 

.....PHI Takes on the VaLUES DPHl,2*0PHI,.,.,LMAXtDPHl, 

AND The L*ST value is Phi » 0 (I.E,. XP(LMAxwl)i-l.), 

.... .similarly G(N,L)rTHE VALUE OF G FOR PHI«L*DPHli 

EXCEPT that G(N,LMAX.1),THE value for PHI.O, 

dimension as 0(NMAXh MAXwO.XPfLNAXAn 

dimension G( 1 0.512), »P (512) »HHN( 10) 
common pi 

COMMON LMaX,LPWR,LMIN,LMaX1,G,XP,X,YS, delta, HHN.NMAX 

common SUPRES.NOcCI.RLLMIN.Ua.LTOTaL 
double precision PI,dPhI,PhI,CNL,CNM1L.CNM2l 
INTEGER TLMaXI .CWTYPE.PP.UA 

logical supres.nocci 

data lHXMX /511/,LMAXR,NHAXR/0,0/ 

If (Lmin.gt,0)go to u 


oo run 


LM!N»1 

U(CWTVPE)a,?,3 
? lFrPP,GT.?)LMrN«PP»2 

IF (M0D(PPf ?) •NE,0,AND,PP.GT,<i)LMINi(PP*i)/2 

GO TO a 
5 UM1N«PP*PP 

a LMAX1 #LMtN+1 

LHAXiMlNOCUMAXitMXMX^ 

IFaMAX.Ett,0)UMAX*LMXMX 

UMAXaHAXO(tMTNaHAX) 

DO 5 K»l,20 
LPmPsK 

TlMAXlaLMAXUL^AXl 
lFfTl«AXUGT*lWAXfl)GO TO 8 

5 LMAXJ.TIHAXI 

6 L^^AXaLMAxai „ 

lFfLMAX,EU,LMAXR ,AND, NMA X , l,£ • NMa XR ) «E TURN 
UMAXRatMAX 
^ NHAXRcNhAX 

or* L 2 «lMAxj /2 

DPHl»Pl/OFtOAT(lMAXn 

.••tfCOMPUTING THE ARRaY ,XP, 

DO lO K»i,L 2 
t3stMAxl *K 

XP{l3)fOCOS(DFlOAT fK)*0PHI) 

XP( K)«*XPa3) 

1ft C.ONTtNUI 
X p a 2 ) s 0 
XPruFtAXl)*! 

...COMPUTING THg array ,G,, aUL THE NEcESSaRY 
’ numbers have been computed and are contained in ,XP,, 
millT IS only necessary to combine them in THE APPROPRIATE 

!i*//.MANNER. 

OO 60 L*l,LMAX 
CNM2 l*i 
NLMUL-! 

INDEXaL 

CNHiL»«XP(INdEX) 



Gn#t)8l^CNH|L 
IF(Nmax,EQ,1) go ^0 feO 
DO SC N»2,NMAX 
NtMj «N*fl 

TNOEXtMODfNLMWLMAxn^i 
CNt»X»>f INDEX) 

I^flNDEX ,EQ, MODfNUHi,Tt.MAXt)fj) CNU**CNL 

6<W^L)«(CNMaL«CNO/2,00 

CNM2u«CNMjl 

CNHiLbCNi, 

SO CONTINUE 

60 continue 

GdfLMAXDaa 
lEfNMAX ,EQ,n CO TO dO 
00 TO N«2 ,NmAx 
GCN ptMAXi )*0 
70 continue 
80 CONTINUE 

Return 

end 

subroutine lNF2YfXpNMAX|HHN) 

,,,THIS subroutine computes the influence functions for YiO, 

pp.THIS version of the program is for the Van SPEIGEL modes, 
DIMENSION MHNfNHAX) 

OOURIE precision PHl,DNpSN,SNMlpSNM2pT0PlpPl 

data topi /,63S6jR772SD0/|PI/S; 1«15R26SJ58R7R5 DO/ 

PHl»DARC05d08lE(^)) 

SNM2«DSIN(PHI) 

HHN(i )*T0PI*fPHNSNM3) 
iFfNMAX ,UT, 2) return 
SNM lgOSIN (2,D0*PHI)/2,D0 
HHN(2)a(PHl*SNMn/Pl 

iFfNMAX ,uT. 5) Return 

00 10 NgJpNMAX 

DNaN 

SNaOSIN(ON*PHI)/DN 

HHN(N)»(SNM2.SN)/Pj 

SNM28SNH1 



SNM^eSN 
in CONTINUE 
Kf TURN 
END 

FUNCTION CEL?(Z»Ai8) 

parameter of The euiptic integrais 

GEDal .-Z 
IFtGFD)l#2#<^ 
y return 

SET result value s OVERFLOW 
? lFfB)i.5,a 
i CEL2»*lfEiP 
return 

a cEL2«1,E3S 
return 

5 CEL?»A 
RE TURN 

COMPUTE integral 
^ GtO«SORT fGEO) 

ARTil , 

AAsA 
AN 9 A 4 8 

NsB 

7 WsW^AA*GE0 
WsW4W 

aakAki 

AARIsARl 

ARUGEO4ARI 

ANaW/ARl+AN 

TEST OF accuracy 
IF f AARI-GEO-1 

fl GEn*SORTfGEO*AARn 
GEQvGEO^GEO 
r.O TO 7 

R CEL2»,TSS3R«2 *AN/ARI 

Return 

Eno 



